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We have calculated the dispersion relation for the high energy cosmic rays (protons, electrons 
and neutrinos) propagating in the Bose gas as well as Bose condensate medium of the pseudo-scalar 
particles. For cosmic rays propagating in the Bose gas, the mass of the particle will decrease and 
this is proportional to the boson density in the medium. If the propagating fermion is massless then 
it will develop an imaginary mass and thus will be absorbed in the medium. But if the condition 
m 2 /m^pa ~ 1 is satisfied, then the cosmic rays will propagate freely without loosing energy in the 
Bose condensate medium. 
PACS number(s): 98.70.-f, 03.75.Fi, 95.35.+d 
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I. INTRODUCTION 



Axions and majorons are the pseudo-scalar Gold- 
stone bosons associated with the Peccei-Quinn and the 
U(1)b-l symmetry respectively and they get mass be- 
cause of the spontaneous breaking of these symmetries 
jj],U. Astrophysical and cosmological calculations allow 
the axion mass range 10 -6 eV to 10~ 3 eV. Axions in this 
mass range would have been produced non-thermally, as 
a Bose condensate i.e., a classical coherent field oscilla- 
tion of the axion field ||. Because the axions are pro- 
duced with small velocity, they behave as cold dark mat- 
ter (CDM), in spite of their small masses. Also the CDM 
scenario for structure formation is at present the most 
promising with axions and neutralinos being the leading 
CDM candidates. If axions provide the bulk of the dark 
matter, they must comprise a significant fraction of the 
dark halo of our own galaxy B. The existence of co- 
herent axion mini clusters has been suggested B. It is 
shown that, axion clumps were formed in the universe, 
when the temperature was about 1 GeV, because of the 
nonlinearity of the axion potential and the inhomogene- 
ity of coherent axion oscillation of the scale beyond the 
horizon ||. 

It is expected that, in the present universe, there ex- 
ist the axion mini clusters and the axion boson stars as 
well as the incoherent axion gas as dark matter comprises 
these coherent axion clumps. Majorons of KeV mass is 
also proposed as the candidate for the dark matter M . As 
the majorons are also pseudo-scalar particle, they can un- 
dergo Bose Einstein Condensation (BEC) and can form 
clusters in the galactic halo. 

The discovery of ultra-high energy cosmic rays with 
energy above Greisen-Zatsepin-Kuzmin (GZK) cut-off 
(~ 5x 10 19 eV) M are particularly mysterious since, 
if the cosmic rays are protons, then photopion produc- 
tion caused by the resonant scattering process with the 
microwave background should result in a rapid loss of en- 
ergy. It is believed that all the high energy cosmic rays 
are AGNs origin. The propagation of these high energy 
particles in the inter galactic medium is also presently a 
puzzling phenomena. 



In this present work we study the propagation of ultra- 
high energy cosmic rays in the bosonic medium and BEC 
medium. The paper is organised as follows: propagation 
of the fermion in the Bose gas is studied in section 2. Sec- 
tion 3 deals with the propagation of fermion in the BEC 
medium. Ultra-high energy cosmic rays propagation, in 
Bose gas and in BEC medium is discussed in section 4. 
In conclusions we briefly summarize our results. 



II. FERMION SELF-ENERGY IN A BOSE GAS 

At finite temperature and density the physical pro- 
cesses take place in a heat bath of particle and anti par- 
ticles and the properties of the test particle deviate from 
its vacuum values. The dispersion relation satisfied by 
the free propagation of the particle is modified, in the 
presence of the heat bath. The dispersion relation of the 
photon and leptons are well studied in the finite temper- 
ature heat bath, particularly in astrophysical as well as 
in cosmological scenario. However to my knowledge the 
particle propagation in Bose gas and in Bose condensed 
matter is not yet well studied. 

The fermion and pseudo-scalar interaction is given by 
the lagrangian density, 



£int = igpi^ibH 1 



(1) 



where 4> is the pseudo-scalar field and ip is the fermion 
field. The self energy for fermion in the presence of a 
heat bath is given by S , 



«£ = 



d 4 k 



{ig p ^iS F {p - k)ig p "f b iD(k)) (2) 



where Sf{p) and D(k) are the fermion and the bo- 
son propagator. In real time formalism of finite tem- 
perature field theory the fcrmionic and the bosonic 
(scalar/pseudo-scalar) propagators are given by 



Sf(p) = {fr + m) 



1 



(p 2 — TO 2 ) 



iT F (p) 



(3) 



and 



D{k)=w ±- T iTB ( k) 



(4) 



(5) 



respectively. The quantities 

T F (p) = 2n5(p 2 - m 2 )n F (p ) 
and 

r B (k) = 2TTS(k 2 ~ml)n B (w). (6) 

Let us consider the propagation of a fermion in the 
medium. The real part of the self-energy term gives the 
propagation of the particle in the medium, and the imagi- 
nary part gives the absorption. Also let us consider that, 
only bosons (pseudo-scalars) are there in the medium. 
Then the real part of the self-energy term will be given 
by 

p v 2 [ d 4 k 5(k 2 -ml)n B (w) 

ReT, = -g p J -jr^-^yp- A + m)- 



(2tt)S 



(p - kf 



(7) 



The real part of the self energy can be expressed as 

-ReE = ap + b jd (8) 

where a and b are two Lorentz-invariant functions and 
u^ = (1,0), pi* = (p°, -p) and k>* = (w, -k). Then the 
coefficients a and b are given by 



[Tp - T u po) 



and 



( ,= lg-i)r.-»T„ 



where T p and T u are given as 



and 



T p = —Tr{£>Reiy) 



Tu = --Tr{/iReZ). 



(9) 



(10) 



(11) 



(12) 



In general the dispersion relation for the fermion is given 
by 

det(^-m-E) = (13) 

and this gives rise to 

(po(l + a) + 6) 2 = m 2 + 2(1 + a) 2 p 2 . (14) 

By substituting eq.(ffl) in eqs.(0) and (O) and the sim- 
plifying them we obtain 



kdk 



T = ^p 

p 32^ 2 |p|7 Wfc 



n B (w k ) [{pi p 2 + to 2 - m 2 ) 



x In 



D 2 + cr 



D ? -a 



2a 



and 
T - 



(2 



ll (kdk 

-vo / o — vP° - Wfe J ns(wfc Jin 



£>2+cr 



Z?2-CT 



where u>fc 



£>2 = Po _ P 2 _ ™ 2 _ "4 ~ 2 PoWfc 



(15) 



(16) 



(17) 



and a = 2|p||k|. For fermions having po and |p| 
>> to, TO^ and |k| the lograthmic term will vanish i.e. 



In 



£>2 + cr 



£>2-cr 



0. 



This gives 



T„ 



9p f k 2 dk 



8tt 2 



U'k 



n B {w k ) 



(18) 



(19) 



and T u ~ 0. Putting T p and T u in eqs.(0) and (|l|) we 
obtain 



-apo- 



(20) 



For exp((3Jk 2 + mV) >> 1 we can write the integrand 
in eq.(]19|) as a sum. Then 



f k 2 dk i ~ r~ , 

a 



^ _ a 2 e -(n+l)x da . 



71=0" 

a ^ i^_i (a(n + 1)) 
~~2 2_^ 



n=0 



71+ 1 



(21) 



where a = to^/3. In the limit a(n+ 1) — > oc the modified 
Besscl function 



K_ x {a{n + l))- X: , 

2a(n + 1) 

Putting this in eq. (|2l|) , we get 

'-2, 



p -(n+l)a 



1 



8n 



(22) 



/ n B (wfc) = > 



-(n+l)a 



>V2/3 4 (n + l) 3 



8o: 

(23) 



Keeping the leading order term in n and for m^P >> 1 
we have 



T P = 



% 



4rri0 



«o, 



(24) 



where 



/ m^ \3/2 



(25) 



is the number density of bosons in the medium. Then 
the dispersion relation in eq.(ll4|) is simplified to 



(26) 



pi - p 2 - to 2 = p 2 (a 2 + 2a). 



Keeping terms up to order a in the above equation the 
dispersion relation is 



pI-v 2 



9 2 P 
2m„ 



n . 



(27) 



Thus in a bosonic medium the propagating particle mass 
square decreases by an amount Sm 2 = g 2 n Q /2m c j i . If 
the propagating fermion is massless then it is Landau 
damped, because it will acquire an imaginary mass. On 
the other hand if the condition m 2 = g 2 no/2m l f > is satis- 
fied then the fermion will propagate like a massless par- 
ticle in the vacuum. 



III. FERMION SELF-ENERGY IN BEC MEDIUM 

For BEC of the pseudo-scalar bosons, the number den- 
sity rip can be expressed as || 



n B (k ) = (27r) J n <$ d (k) 



(28) 



where no is the particle density in the condensate 
medium, with zero momentum. Then putting eq. (|2§|) 
in eq.(0) we get 



T n = 



Oo - P - Pom 4 



n 



P ~ Up ( n 2 _ „2 



(Po - P 2 + TO I ~ 2p O "i0 - m 2 ) 2m 4 



(29) 



and 



Tu = -3p" 



(po - w) 



»o 



(Po - P 2 + mi - 2pom<p ~ ™ 2 ) 2m ? 



(30) 



Then putting the values of T p and T u in eq.(|9|) and (|l 
we obtain for the BEC medium 



C« 



(po — m^) 2 — p 2 m 2 



and 



b = - 



am^,, 



(31) 



(32) 



where C p = g 2 no/2m t f,. For BEC medium b is propor- 
tional to the pseudo-scalar mass whereas for Bose gas it 
is proportional to the energy of the propagating fermion 
as shown in eq. (p0|). Putting the value of b we obtain 



2„2 



(po(l + a) - arritf,) =m + (1 + a) p 



(33) 



Keeping terms up to order a the dispersion relation will 
be 



9 

Pa 



1 2 

p — m 



(pi 



2C p m ( j,po — 2m 2 C p 

p 2 - m 2 ) + to 2 - 2p Q m ( j, ' 



Writing p 2 , — p 2 — m 2 = Sm 2 we get 
2C p (m ( / ! po - m 2 ) 



Sm 



5m 2 + (to 2 , - 2p m tt> ) ' 



(34) 



(35) 



It is a quadratic equation in <5to 2 , so 6m 2 has two solu- 
tions. But the physical solution is the one, which goes to 
zero for C p goes to zero. So the value of Sm 2 is given by 



Sm = 



2p - TO0 



2C p (po 



For pq » rriif, the above relation will be 



Sm 2 



Pa 



Pa 



TO TO TO^ 

TO0 2p 2 



(36) 



(37) 



Now the dispersion relation for the fermion propagating 
in the BEC medium is 



Pi 



9pno 
2mA, 



+ 



m,pp 2p\ 2p J 



2+) (38) 



where is uq is the number density of the pseudo-scalar 
particles in the medium. Comparison of cq.(p7|) with 
eq.(|3S|) shows that, we do not have a pure BEC medium, 
because the second term in the right hand side of eq. (|3q) 
i.e. g 2 no/2m,p comes from the Bose gas contribution. 
Again the eq.(p8|) clearly shows that, for a massless 
fermion propagating in the BEC medium, the particle 
will develop a negative mass square value of magni- 
tude g 2 no(l + m ( p/2po)/2m c / > and will be absorbed in the 
medium. Now let us compare the contributions of differ- 
ent terms (within the bracket of eq.(p8|)). For high energy 
particle po >> m^- So the term m ( / ) /2po almost does not 
contribute. For particle with po >> to the term m 2 /2po 
is also small. So only the first and the second terms with 
in the bracket are contributing. Thus we have 



Pi 



gp"0 

2ms 



to 2 \ 



(39) 



Let us consider different situations, which might arise 
depending on the mass and energy of the propagating 
particle. For massive fermions if to 2 /to,/, >> po condi- 
tion is satisfied then the particle mass will be increased 
by an amount g 2 n^m 2 j2m\p§ and the BEC contribution 
will dominate over the Bose gas one. On the contrary if 
to 2 /to0 << po then, the particle mass will be reduced 
by an amount g 2 no/2mA, which is purely due to the Bose 
gas nature of the medium and BEC contribution is neg- 
ligible. For po — m 2 /m^,, the particle will propagate like 
a free particle in the BEC medium without loosing any 
energy. 



IV. COSMIC RAYS PROPAGATION 

Active galactic nuclei are found to be the origin of the 
highest energy cosmic rays, like photons, protons, elec- 
trons and neutrinos p0[ . The recent discovery of cosmic 
rays events above Greisen-Zatscpin-Kuzmin (GZK) cut- 
off (~ 5 X 10 19 eV) § by the AGASA Q, Fly's Eye ||, 
Haverah Park |L|] and Yakutsk collaborations |L4|] is an 
outstanding puzzle in cosmic rays physics. If these cos- 
mic rays are protons with energy more than GZK cut-off 
then, interaction with the cosmic microwave background 
radiation (CMBR), would cause rapid loss of energy by 
photopion production and consequently depletion of the 
observed flux of these particles [J7| . For every mean free 
path ~ 6 Mpc of travel, the proton loses 20% of its energy 
on average Q . Since AGNs are hundreds of megaparsecs 
away, the energy requirement for a proton which arrives 
at earth with a super GZK energy is extremely high. 
Presently we do not have a reliable theoretical model to 
explain the origin and propagation of such high energy 
cosmic rays in the galactic medium. On the other hand 
ultra-high energy neutrinos will escape first and bring us 
first-hand information regarding the source. Also these 
AGN neutrinos can be used to study the properties of 



neutrinos themselves |10 



Here we will see how the propagation of these ultra- 
high energy cosmic rays (protons, electrons and neutri- 
nos) are affected by axion and majoron dark matter in 
the galactic halos. 



A. Bose gas 

Let us consider ultra-high energy cosmic rays proton 
propagation in the axion gas. Then for proton-axion 



coupling g p ~ 10 



m$ 



10~ 3 eV § and 



axion number density uq ~ 3 x 10 13 /cto 3 we obtain 
g 2 no/2m a ~ 1.2 x 10~ 18 eV 2 and this is much smaller 
than proton mass square mi. Also this term is much 
smaller than electron mass square m 2 (we assume the 
electron-axion coupling same as the proton-axion cou- 
pling) . So the axion gas has no effect on the propagation 
of the cosmic rays protons and/or electrons. Secondly 
let us consider the neutrino propagation in the majoron 
gas, then taking neutrino-majoron coupling g p ~ 10~ 4 , 
rnj ~ 1 KeV [Q and the majoron density same as the 
axion density, we obtain g 2 riQ/2mj ~ 1.2 x 10~ 12 eV 2 . 
For non-zero neutrino mass the majoron gas has also no 
effect on the neutrino propagation. 



B. BEC medium 

For proton propagating in the axion Bose condensate, 
we have m 2 /m a ~ 10 21 eV which is above the GZK cut- 
off. For proton energy p << 10 21 eV the proton mass 



will increase by an amount g 2 nom 2 /(2pom 2 ). Similarly 
for electron propagating in the BFC of axion, we have 
to 2 /to0 ~ .25 x 10 15 eV. In this case also for electron 
energy po < 0.25 x 10 15 eV the electron mass will in- 
crease in the medium and for po >> 0.25 x 10 15 eV the 
mass will decrease. If m 2 /m a po — 1, then the Bose gas 
contribution will be cancelled by the BEC contribution 
and the proton will propagate like a free particle in the 
vacuum, without loosing its energy. 

For neutrino propagating in the majoron condensate 
medium, we have m 2 /mj ~ 10 _1 eV, for taking neu- 
trino mass m„ = 10 eV and raj — 10 3 eV. The cos- 
mic ray neutrinos have energy as high as 10 21 eV. So 
m 2 /mjpo — 10 -22 ; hence only the majoron Bose gas will 
contribute for the neutrino mass and the neutrino mass 
term in eq.(|38|) will be (m 2 — g 2 n^/2raj). This reduction 
in neutrino mass is very small to account for. 



V. CONCLUSIONS 

We have calculated the dispersion relation for the 
fermion propagating in the Bose gas and the BEC of the 
pseudo-scalar particles. It shows that, if the fermion is 
massless then, it will develop an imaginary mass both in 
the Bose gas and in BEC medium, which is proportional 
to the number density of the bosons (Bose condensate) 
in the medium. For cosmic rays of energy below ra 2 /ra^, 
the mass of the fermion will increase, on the other hand 
for energy above ra 2 /m^ the mass will decrease. Inter- 
esting situation arises, when m 2 /ra^po ~ 1. In this situ- 
ation the Bose gas contribution will be cancelled by the 
BEC contribution and the particle will propagate freely 
without loosing its energy. But we found that the contri- 
bution of the axion gas or BEC of axion to proton mass 
is very small. Similarly majoron has also very small con- 
tribution to the neutrino mass. This small contributions 
are solely due to the low number density of these particles 
in the galactic halo. 

It is a great pleasure to thank Dr. S. Mohanty for 
many useful discussions. 
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